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ABSTRACT 


The hydrodynamic stability of an underwater station was studied 
using wind tunnel models. Rigid model force and moment coefficients 
were obtained for various orientations of the station in three 
dimensions. The unsteady motion of the station in cross currents 
was investigated with a dynamically similar model. Color 8 mm 
pictures and amplitude-time traces of the unsteady motion were 


obtained. 
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LI. INTRODUCTION 


The purpose of this study was to investigate the hydrovynamic 
stability of a manned underwater station when the station is fully 
submerged and subject to cross currents. The concept of a manned 
underwater station, or MUS, has been investigated by General Dynamics 
Corporation, Electric Boat Division, under contract for the U. 5. 
Naval Civil Engineering Laboratory, Port Hueneme, California (Ref. 1). 
The present study, which is in support of the earlier MUS work, was 
sponsored by NCEL under work order WR-8-0128 and was monitored by 
Dr. Cheng Lung Liu of NCEL Ocean Engineering Department. 

The object of the MUS project as stated in Ref. 1 is to develop 
a deep ocean submergence habitat which will accommodate a five man 
Eo for an operational period of thirty days. The station is to 
be capable of submergence to depths of 6,000 feet. 

As envisioned by General Dynamics, the station would consist of 
two circular cylinders of unequal diameters and heights. The smaller 
diameter cylinder is capped with hemispherical ends and will contain 
a nuclear reactor. The larger diameter cylinder is capped with 
spherical ends which are fitted to the cylinder by conical, tangential 
shells. (Fig. 1) The habitat cylinder will have five decks and is 
the working space for the crew. 

Since these cylinders are designed only for a vertical orientation, 
a large ballast tank completely encircles both cylinders. The design 


of this tank is based in part on weight and balance considerations. 


1s 


The station is not equipped with a propulsion device although 
it does have ballast and buoyancy control equipment. The concept 
of the station is somewhat unique and the MUS may best be described 
as an elevator. In concept, the station is to be towed on the surface 
to a pre-selected site by a Support vessel. When the operations 
site is reached, the Oticn will be manned and an emplacement 
procedure will commence. The MUS will utilize an anchor and cable 
arrangement to pull itself to the bottom in bootstrap fashion. It 
is equipped with a winch (located on the main ballast tank) and 
attendant systems for control. 

When submerged, the MUS will maintain a slight positive 
buoyancy at all times. The anchor will first be lowered to the 
ocean bottom while the MUS holds on the surface. Once the anchor 
is in place, the buoyancy of the MUS is reduced and winchdown begins. 
As the winch retrieves the cable, the MUS will haul itself to the 
bottom. The recovery cycle, i.e., the ascent to the surface, is 
just the reverse of this procedure, 

It is primarily with the emplacement and recovery operations that 
this study is concerned, When the MUS is in the mid-water column, 
it may well be subjected to a cross flow caused by deep ocean 
currents. Evidence exists to expect such currents throughout the 
water column, having a magnitude of 0,1 to 0.5 knots. 

If such currents were to exist, it would be reasonable to expect 
the MUS to be forced into motion until equilibrium is attained. 
Such motion could be translational or rotational and perhaps 
oscillatory in nature, Since the MUS is restrained only by a single 


cable (to the anchor), the only other stabilizing elements are: 
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1. The buoyancy couple caused by the displacement of the meta 

center relative to the center of gravity. 

2. Hydrodynamic stabiliaation, i.e., the geometry of the station, 

including such extremities as the legs and pads. 

3. The inertial resistance to motion (moment of inertia) 

Another aspect of the MUS recovery cycle which necessitated a 
stability investigation was that of emergency ascent. It was proposed 
that, in the event of an emergency such as cable jamming or bottom 
entrapment of the support pads, a free ascent could be accomplished 
by cutting the cable (with explosive cutters) and/or jettisoning the 
support pads and legs. Since these are in fact stabilizing devices, 
such an ascent could only be feasible if the resulting configuration 
of the MUS provided sufficient stabilization. 

The objectives of the experiments reported herein were twofold, 
namely: 

a. To establish by means of rigid model wind-tunnel tests the 

pertinent dimensionless force and moment coefficients for 
the MUS at various orientations, and 

b. To determine the dynamic behavior using a "partial mode" 

dynamic model, and in particular to note the existence of 
any self-induced and/or vortex-induced oscillations. 

In this investigation, the MUS was represented by an 0.0209 
scale model. Reynolds numbers were matched between model and prototype 
so that aerodynamic measurements on the scale model were applicable 
to the hydrodynamic traits of the prototype. 

One should note that the MUS configuration is still in the early 
design stages, hence the results obtained in this study should have 


as‘a goal the possible application to other configurations. However, 


li, 


accurate forecasting tends to be uncertain because the bluff-body 
features of these configurations do not interact in a linear manner. 
The testing accomplished in this study was relatively easy, and other 


configurations could be studied using a similar approach. 
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II DISCUSSION 


A. DESCRIPTION OF MODELS AND EXPERIMENTAL PROCEDURE 

The investigations described herein upon the 0.0209 scale MUS 
models were performed at the Naval Postgraduate School (NPS) low- 
speed wind tunnel facilities. The force and moment measurements 
were conducted on the "rigid" model in the Aerolab 32 in. x 45 in. 
wind tunnel while the scaled dynamic model was tested in the West 
Coast Research 48 in. x 60 in. wind tunnel. 

The "rigid" model was linearly scaled and was fabricated from an 
overlay of Honduras mahoganyy upon an aluminum spar or reinforcing 
plate (Fig. 6). The purpose of the spar plate was to make possible 
the attachment of the support arms and tail sting at prelocated ox 
orientations of the model. It was possible to position the model 
in a multitude of ways upon the conventional three strut support 
system so that representative yaw and pitch attitudes could be 
obtained for an angle range of Siee and 90 respectively. Figure 
5 presents a sketch of the model including pertinent dimensional 
data relative to the force model, while Figures 11 to 44 show the 
model installed in the wind tunnel at various orientations. The 
model was tested at dynamic pressures from 3.85 to 57.8 psf (in air) 
which correspond to prototype velocities in sea water from 0.08 to 
0.30 knots (refer to Figure 52 for conversion from model to prototype 
operating conditions), 

The forces and moments were measured upon an Aerolab three- 
component beam balance (Fig. 10) which had the capability of measuring 


lift, drag, and pitching moment relative to the wind axes. 


i 


The bollancewibimicts were: 


Component Limits Sensitivity 
kee roves Somlibe = cOMOmOmmne 
+ + 


Drag force 50 lbs 0.005 lbs 


1+ 
Pe 


Pitch moment 1000 in-lbs 0.050 in-lbs 

Care should be exercised in interpreting the meaning of lift and 
drag forces and pitching moments at some of the extreme model positions, 
such as when the model has been rolled 90 degrees on the balance, 
since in this attitude pitch angle changes actually corresponded to 
model yaw angle changes (Figs. 43 and 44). In this extreme case, 
lift force and pitching moment correspond to side force and yawing 
moment respectively. It is unfortunate that the balance system was 
not readily usable as a six-component type since then it would have 
been possible to convert the wind axes results to stability or body 
axes as is frequently done in aircraft analysis. 

Data corrections were made for weight tares, while aerodynamic 
tares to account for the influence of the support system were 
approximately estimated by measuring the balance force and moments 
without a model present. Since the model was a relatively bluff body, 
as compared to a streamlined aircraft configuration, these 
approximations were considered as being reasonable. 

The dynamic model tests may be described as "partial mode" tests 
since only selected degrees of freedom were investigated. The modes 
investigated included lateral and axial translation and yaw rotation, 


The tests were performed upon two 0.0209 scale models, one made from 


balsa and the other from sugar pine (Figs. 7, 8, and 9). 
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The balsa model had as its purpose an initial observation of the 
model's dynamic behavior and an experimental confirmation of the 
yawing moment static trim points. The full scale weight and weight 
distribution table for the MUS from Ref. 1 was used as a basis for 
design of the pine model. By using a component approach, the full 
scale weights and weight distribution were scaled to obtain the desired 
mass moment of inertia for the pine model. Removable lead inserts 
were utilized to allow a variation of mass and mass moment of inertia. 
The balance of the model about the desired center of gravity location 
was checked experimentally, while the mass moment of inertia was 
determined by utilizing the principle of the torsional pendulum. 

An object having a known moment of inertia, namely a solid steel 
cylinder of given weight and geometry, was suspended on an 0.037 

inch diameter steel piano wire of known length. The oscillation 
frequency was then measured which in turn allowed the determination 

of the wire's shear modulus of rigidity. From the physical properties 
of the wire (diameter, length, and shear modulus), and the measured 
free vibration oscillation frequencies, it was possible to determine 
the mass moment of inertia of the model. Table II lists the physical 
properties of the dynamic models tested. 

A narrow diameter tube was inserted along the vertical axis 
(through the assumed prototype center of gravity) in order to shield 
the dynamic model from the support wire. The model was suspended on 
an 0.037 inch diameter steel piano wire which was fixed to supports 
external to the tunnel. (Figs. 56 through 59) The wire was preloaded 
to 50 lbs. tension by means of turnbuckles, and the tension was 
frequently checked by using a standard rigging tensiometer. Initially, 


only the taut wire support was used, but later tests were conducted 


“zd 


with the translational degrees of freedom removed by means of local 
tie-down constraints. In addition, the last sequence of dynamic 
model tests were provided with an electrical angle sensing potentiometer 
in order to allow time-history traces to be made of the angular motion. 
(Figs. 60 and 61) 

Motion picture photography of the dynamic model tests was accom- 
plished by using a mirror arrangement. An adjustable mirror was 
mounted to the tunnel wall above the model so that overhead photography 
was possible. The later tests involved simultaneous photography 


and electrical data collecting. (Table VI) 


B. SCALING CONSIDERATIONS 

The scale of the models (0.0209) was established on the basis of 
the largest model compatible with the wind tunnel test section size 
and velocity range in combination with a Reynolds number match 
between model and prototype. 

For the prototype design with an ocean current of 0.3 knots in 


salt water, an estimate of full scale Reynolds number would be 


hoes 
Noe = Vedp = (0.3 KTS x 164 see unas. St) 


ase 


sai ih inn 7: 


so 
Nee = 4-45 x /O (prototype in salt water) 


The above Reynolds number is based upon a mean diameter of 14,93 
feet. An 0.0209 scale model (A = 0.0209 = scale ratio) would satisfy 


the prototype Reynolds number condition above at a tunnel velocity of 


—4. oe = 
i Um Nee = L,S68 xX /O fee (4.45 x 10 ) = Aeon. Ts 
dm ( Om@inmetia 


ap 


which corresponds to a tunnel dynamic pressure of 60.3 psf. The ratio 
of model to prototype velocities is 443 under the assumption of 
equivalent Reynolds numbers. 

The mass moment of inertia scaling of the dynamic model was 
established by the free vibration equation in the yaw mode. As 
an approximation, one may consider the slope of the Cy vs Y curve 
to be linear in the immediate neighborhood of a zero moment static 
trim position. 
The total yawing moment may be approximated by the sum of the 


linearized static term and a term due to the dimensionless angle 


ate (ES. 
, 
Ps eee lov, | ¥ cn) 


Multiplying Cy, by (gad ) to obtain physical units, then the free 


vibration equation in yaw is: 


Ts Ly - (gaa) (& ous, SY ~ (gRe) Cw, Y=0 


re) 
where: i = Mass moment of inertia, taken about the # axis. 
Dividing through by T5,.and introducing the dimensionless time, t 
where: f= a 
(4 
ZV 
then a 
(20) d%, ~ BAL Cw, ci. o m, bbe y=o 
GS Le Te dt i 
ez ies 1CD 


The enforcement of Reynolds number equivalence insures that the 
dimensionless yawing moment coefficient Cv is ident teaver 


both model and prototype. In order to achieve full dynamic similarity, 


Zz 


the nondimensional, natural frequency of the model must be identical 
with that of the prototype. This argument then defines the scaling 


law for the mass moment of inertia as shown below. 


3 c 
sry mes pad (¢c St la Pad Ca = 
d he a Ny) its 6 a ) : eee 


This may be compared to the familiar vibration equation 





las +26 W Sot = (o) = 
A a= EM lal — 


In order for the equation of motion to describe concurrently the 
behavior of both model and prototype, each term must equate in 
model to prototype conversion, therefore 


S, : ¥ mode - Len a 


Prototype 


or 3 3 
I nee (ow, Yin = Pp et F Cn yy, VP 


YY\ 


since Cia en a Cny 
ers = Ve 





3 
Ae d*p 
then 
5 
aE 
Pres Tee, 
Os Ay 





em = Fea a 
- 
= 6 oie come 


if the inertial scaling of above were used, then the non-dimensional 


frequency would be alike between model and prototype, i.e., 


24 


Similarly, a check of the damping terms: 
3 
24 — ] . ie are 


will yield the same result by the same argument. 


Unfortunately, such scaling was not physically realizable as it 
would demand a prohibitively light construction of the model. 
Returning to the non-dimensional natural frequency, which may 


be identified in eqn. (3b); 


3 
— 2 
WW, = 7 . PAS Cy 
tas ¥ weseee (4a) 
— 
nm 2 (YN (Te 
(Coy, > Fe np 
if the mass moments of inertia were set as: 
Titans, 2 
Ip 7 


then the dimensionless frequencies won't correspond, but will scale 


as: — 2 - 
ok 7 


a 


(W," Dp tp 
Considering that dimensionless frequency is by definition: 
Wn, = Zy Yn 
one readily obtains that: 


Yr, = (SY(EV(FE (on, 


Pa 2) 


ah, 


After substitution of the scaling parameters (cf. Table I), the true 


frequency scaling becomes: 


(np = 0.00134 Wn), 


For example, if (ae =e 70 EPS , then 


Se cycles 
(Wn>d— = 0.00134 (070 <Ps)(60 mn) = 0.0562 iin 





and the prototype period would be 18.1 min/cycle. 
Using a quasi-steady argument, enforced by the case where the 
damping ratio e is zero, one may use the above scaling law 
to convert model frequencies to prototype frequencies for the yaw mode. 
In a similar fashion to the above analysis concerning the yaw 
mode, one may consider the translational mode free vibration equation. 
In the lateral translational mode, a similar analysis to the yaw 
mode scaling may be made starting from the free vibration equation, 
mY - "ea (BAD 4 + 4 aC) a 
The damping derivative, Ce. » may be approximated by the quasi- 
steady assumption in a manner similar to that used by Parkinson, 
(Ref. 2) when he considered the limit-cycle solution for the galloping 
transmission line problem. Using the sign convention of Figure 2, 
one obtains that: 


C = + 


‘aide 


Y 
since v= nF YK, 


Substitution of this approximation into equation (7) yields that: 


siafilers'« COD) 
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It iS convenient to convert equation (8) to dimensionless time, L ; 
av) vy 

and dimensionless lateral translation, 4 , where a = Y : 

Direct substitution plus a slight amount of algebraic simplification 


yields: 


d*G _ ¢Ad 


omen: ae nae 
ag 4m ry Te + E ¢ = 


a O 

= 
4v a; «rece 
Equation (9) may be compared with the second order differential 
equation; equation (4), in order to identify that in the lateral 


translational mode, the non-dimensional frequency is: 


aR, = (4) 


where ae = ( % 


For the model installation, the spring constant, k, is proportional 
to the tension ih in taut wire support as may be: seen in the rear 


view sketch below: 


— or GR N SN NNN 
) Nevtra| deSlected 


laterally 








An approximation for the tethered prototype vehicle may be made on 
the assumption that the cable tension is equivalent to the buoyancy 


force. For this case, one would obtain 


Soe 
Ro Lp 


2/ 


where the tether cable length, S » 1s an equivalent value which 


may be estimated. 
The initial consideration will be to find the true frequency value 


of the prototype on the assumption of equivalent non-dimensional 


frequencies, followed by a correction for differences in cable tension. 


Cena 


il 


F the se of — 
Ons ca O (2), 


if 


it follows that: 
(ae 


(? 


(Se . (= es in 


= | Q 
(Un ¢ 7\450 ( nen 





Since 


Next, we shall investigate the implications of cable tension. 


Bch = tes , tt Beollows thaw 
4 = 
(alee Sealy oY ak “age 
+ Ie; | ll 
(Bes 2 Gey Cae) > 
kw, ¥ 1. hk 
i 4600 ua 





Ae = 600 ft., 


cable lenghbs are approximately given by 


ie 
cable tension is 50 lbs., then 


h= 6 ft, and model 
Tipeinee 07 enema arias 
Qo 2300 Sevw 
or 
S 2.18 |b 
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Actually it is envisioned that the prototype buoyancy force will be 


on the order of 100 lbs., hence the corrected frequency would be 


ah 
mS 100 2 
er = ee (2), 


corr 


cep “Celie 
(Qa), a 3170 Cas 


Therefore if the lateral translational frequency of the model were 
about 1 cps, then the representative prototype period would be 
about 53 minutes per cycle. 

In retrospect, the dynamic scaling problem would require 
simultaneous scaling of mass, mass moment of inertia about the 
yaw axis, and cable tension. Since this procedure was not initially 
followed, estimates were made to account for the relative effects 


of the mismatches, 
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A, RIGID MODEL DATA 

Various rigid model plots are included (Figs. 45 through 55). 
Test runs 1 through 37 indicated that the hydrodynamic coefficients 
are essentially constant throughout the Reynolds number range under 
consideration. A large number of Coefficient-Reynolds number plots 
were obtained (see run schedule, Table VI), and representative 
plots are shown in Figures 46 through 48 inclusive. The largest 
values of the hydrodynamic coefficients that were measured, with the 
corresponding full scale forces and moments, are plotted for 
comparative purposes in Figure 52. These forces, as compared to 
the buoyancy force of 100 lbs., are considered to be significant. 

A drag coefficient contour plot is included in Figure 45. The 
drag coefficient values along the ordinate and absisca of this plot 
were measured in five-degree increments (See also Figs. 49 and 54), 
while other values were measured in combinations of thirty-degree 
increments, 

The plot of yawing moment coefficient versus yaw angle (Fig. 53) 
shows that the MUS has two stable trim points in yaw, at zero degrees 
and 180 degrees, while two unstable trim points exist at 90 and 270 
degrees of yaw. It should be emphasized that these are zero yawing 
moment positions, while the slope indicates the degree of stability. 
Also note that stable trim did not occur at the orientation for 
minimum drag. 

Direct measurement of the yawing moment vs yaw angle has 


application to defining the aerodynamic restoring moment term in the 


30 


yaw-degree-of-freedom dynamic equation. However, means do not exist 
for direct measurement of the damping term due to yaw rate. The 
damping term must be deduced indirectly by observing the time history 
of decaying yaw oscillation. 

The plot of side force versus yaw angle at zero pitch angle 
is shown on Figure 55. An approximate analytic expression for the 
variation of side force coefficient is: 

Cp 0.38 Sth ¥ 
However, the actual plot shows many irregularities in the curve 
which may be related to the complex flow interactions between the 
various bluff-body components that comprise the MUS. Especially 
noteworthy is the rapid unstable trim points, ¥ = too . 

It should also be noted that the above expression for side force 
coefficient will result in the damping term of the lateral trans- 
lational degree of freedom having a negative value. The consequence 
of this feature was noted in the dynamic model as a limit cycle 
type of motion. 

B. DYNAMIC MODEL RESULTS 

A color 8 mm motion picture film:was prepared as a supplement 
to this report. Part one of the film shows the experiments with 
the single wire arrangement, and careful observation will show the 
lateral motion beating phenomenon. In this configuration the model 
had degrees of freedom both in yaw and lateral translation. The 
combination of an apparently lightly damped yaw mode plus a 
negatively damped lateral translation mode produced a complex 


response with an apparent indication of beating due to the closeness 
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of the two frequencies. The lateral excitation caused the self- 
induced nature of the oscillation in yaw, mae this did not appear 
when the lateral constraint was imposed. 

Part two of the film shows the instrumented and constrained 
model experiments. No lateral motion occurred and the damping 
ratio was small but finite until a dynamic pressure of five or six 
psf was reached. Strip recorder charts are included. (Fig. 67) 
Dynamic model data are summarized below in Figs. 62 through 66. 
The predicted frequencies for the model, shown as estimates on 
Figs. 62 to 64, are based upon the recognition of frequency from 


equations (2) and (4a) as: | 
2 
_ 2v— _ (GA sc 
U4, la AS ji Te, Nw ) 


Although the predicted frequencies are about 10 to 15 percent higher 
than measured values, the trends with dynamic pressure are proper 
and add confirmation to the concept of an equivalent hydrodynamic 
torsional spring term, It is quite possible that the actual observed 


frequencies reflect rate dependence of the hydrodynamic torsional 


spring term on the bluff body. 
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IV. ANALYSIS 


A. SPamec STAG Derry 

The MUS stability in the pitch mode may be interpreted from Figure 
50; i.e., the pitching moment coefficient versus pitch angle plot. 
This plot was obtained only for zero degrees of yaw. At minus ninety 
degrees pitch angle, i.e., the ascending mode, the MUS will trim 
but be statically unstable. While the pitching moment itself is 
zero at OX = -90., any increase in pitch angle results in a positive 
pitching moment which tends to increase the pitch angle further. 
Similarly, a decrease in pitch angle results in a negative pitching 
moment. At plus ninety degrees pitch angle i.e. the descending mode, 
the MUS is again statically unstable. In the neighborhood of zero 
degrees pitch angle, the MUS has a positive pitching moment coef- 
ficient and hence is not trimmed. Therfore, in the absence of the 
cable constraint, the MUS is statically unstable in the pitch mode. 
The balsa wood dynamic model was suspended in a manner so as to allow 
a pitch degree of freedom for one set of exploratory investigations. 
Its dynamic behavior tended to confirm the preceding conjectures since 
the stable trim points appeared to be weak and the model was quite 
prone to rotate in pitch. 

In the vertical translation mode the MUS is statically unstable. 
Referring to Figure 51( x =0, Y — 0, SVS) , in a direct cross current, 
the MUS has a negative lift, which will tend to cause the MUS to 
descend. This in turn would increase the effective pitch angle, 


and result in a more negative lift. 
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Figure 53 shows that the station is statically stable in the yaw 
mode (at zero pitch) at yaw angles of zero and 180 degrees and 
statically unstable in yaw at 90 and 2/70 degrees of yaw. 

Figure 55 shows that the station is statically unstable in lateral 
translation at yaw angles of zero and 180 degrees. If a lateral 
translational velocity occurs, an effective yaw angle is introduced 
in the nature of sideslip. When this occurs, the side force is 
such that the lateral translation tends to be further increased 


(forced). 


B, DYNAMIC STABILITY 

In the absence of constraints, the experiments indicated that the 
MUS would be dynamically unstable in all translation modes and in 
the pitch mode. 

The yaw mode appeared to be dynamically stable, as the oscillation 
amplitude either decreased in time or remained constant. Negative 
damping ratios were not observed for the range of tunnel speeds used, 
but it is possible that higher dynamic pressures could have induced 
negative damping. 

The scaling relations for the yaw mode and the lateral translation 
mode indicate that the natural frequencies of these modes are very 
widely separated at full scale. This would indicate that the coupling 
tendency observed in the model experiments would not occur for the 
full scale MUS. Thes most important conclusion from scaling 
considerations is that an oscillation on the model in a wind tunnel 
on the order of a second per cycle period would scale to the prototype 


vehicle in water with a period on the order of 20 to 50 minutes per 
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cycle. Therefore it is quite possible that unusual oscillational 
behaviors in the tunnel would not be perceptible to the crew enclosed 


in the MUS because of the low velocities involved. 


C. COMMENTS ON MUS CONSTRAINTS 

The experimental results indicated that the proposed free ascent 
mode would be practicable. Considering an assumed stabilized condition 
where the MUS would be floating freely in any existing current, with 
cable and pads jettisoned, and in the absence of an unequilibrated 
buoyancy force, no hydrodynamic force would be generated. If the 
positive buoyancy force of approximately 100 lbs. were considered, 
the MUS would be forced into vertical translation and hydrodynamic 
forces would be generated. In this orientation, the wind tunnel 
axis position of minus ninety degrees of pitch angle would be 
applicable (Fig. 2). Assuming that a constant side current existed 
throughout the water column, the MUS would begin to rise, sensing 
only the hydrodynamic forces corresponding with Of = 490° 3 
However, in this orientation, both a drag force and a lift force could 
be generated. If one assumes that the MUS could achieve a vertical 
velocity of 0.2 knots, from the wind tunnel results, the resulting 
drag force would be 100 lbs., and the lift force would be zero 
(Figs. 49, 51, 52). Actually, the resulting drag forces would limit 
vertical ascent speed to 0.2 knots, if the 100 lb. buoyancy force were 
maintained. This speed limitation may be verified by cross referencing 
Figs. 49 and 52. The drag coefficient of 0.7 will result in a drag 
force of 100 lbs. at 0.2 knots, therefore the drag force will 
equilibrate the net buoyancy force at 0.2 knots. A check of the 


pitching moment vs. pitch angle plot (Fig. 50) for this orientation 
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shows that the MUS would be at an unstable trim point, that is the 
pitching moment would be zero but the variation of moment with flow 
angle would not produce a restoring moment. As long as the flow 
vector remains constant at od = a5 wwthe «MUS ala remain upright. 
If the flow vector is disturbed from & = -90° er monemiikely, if 
the MUS begins ascending from an off-vertical orientation, the 
resulting pitching moment will be in a direction tending to over- 
turn the MUS. This pitching moment would be on the order of 300 
ft-lbs. maximum (for a 0.2 bee eae speed and a 20° off-vertical 
orientation). Considering the design (buoyancy couple) stability, 
this slight pitching moment should present no difficulties (ref. 1). 
That is,@the loea@tiongot the center of buoyancy above themecmier 
of gravity results in a stabilizing moment. 

Considering jeheweases where the MUS,is tetheredsbypthe icable am 
the mid-water column, if side currents exist with a speed of 0.3 
knots, a pitching moment would be generated on the order of 1300 
ft-lbs. maximum. This would place a side load on the cable attach- 
ment point, as, the cable would equilibrate the pitching moment, In 
addition, a yawing moment would exist, having the same order of 
magnitude. The yawing moment would apply a torsional load to the 
cable or to the cable attachment point, depending on the actual 
design, however, the vehicle would tend to align itself to the 
Stable trim point corresponding to p= Olor 180. 

The oscillatory behavior of the model appears to be of no 
practical consequence to the full scale MUS and probably would not be 
noticable by the crew. Considering the proposed emplacement procedure, 


since thegprincipal static .equulmbrabino Torce that existsijto. counteract 
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drae forces (im side currentsjmis the 100 lbs. buoyancy teree, ene 

wind tunnel measurements indicate that a forced descent would result 
with current speeds as small as 0.2 knots. At this speed, the resulting 
drag force would be equal to the buoyancy force, and therefore, static 
equilibration would not be possible. 

The restoring force F. is equal to: 


F =T sing 


x 

Where T = buoyancy force = cable tension 

and y\ = deflection angle of the cable 
from the vertical 

af T = 100 lbs. 

for a current speed of 0.15 knots and Cy = 0.7 


D = 60 Wes, (Figy 52) 
Sim = 106 
yy = 37° 
or in other terms, a horizontal translation of 3,600 ft. and a 1,200 
ft. descent would occur before static equilibration if 6,000 ft. of 
cable separated the MUS and its anchor, For this reason, it is 
recommended that a higher buoyancy force be considered. 

The foregoing discussion on forced descent is admittedly overly 
simplified and does not take into consideration the special effect 
of catenary behavior of the cable itself. This topic would merit 
individual study and analysis before a complete discussion of 
forced descent could be made. It is also recognized that increasing 
the buoyancy force will create a serious problem involving free 
ascent terminal velocity. In this respect, a compromise would be 


involved between stability considerations:and ascent problems. 
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TABLE I 


VALUES OF SOME PHYSICAL CONSTANTS 


GnaracteelsuLre area (Model yew mete cere «6 @ ee eS ee 0.546 fe 
Chearaccerisitemarea “(prowemy pe mmc ws ss ae se 6 8 1248 35e 
HETeCiit On wid@clcimmts” 900s s a ahr Te ce 6 elicits ce” » slime 1. Oe ae 
He lo htmGieED EOEOLY PC aes s <5 em ewe ill’s” vol Mellie et cl Uemilemncecnn. 4.7 Om te 
Characteristic length '(medel)y Tere. . . es be es 367 Dal 
Characteristic eet. CDIEOEOLY DE) aceaemeemr es « Cemmeh ve iemierts 14.93) fe 
Scale-rarioeare. He. . OR ore. 0.0209 
a ea a Ry 
ea eee cm. kk. 9.0 x 10° 
Pie. . SRP! nox 4,0. sand 
me, saltwater eee ee ee A. ek ee SO 107* fraveee 


Vv, ain e e e ¢ e e e e e e ¢ e « * ® é e ® e e 1.53 x 1074 feaneee 
& 


Pf : 
Pp ° EGGS Bemeek . «oes suo « comme 4 ho 23 xe 


Dimensionless velocity,V. ee ee ee eee 443 
* 
Note: The characteristic length was arbitrarily chosen to be the 


average or mean of the diameters of the five principal 


components which comprise the MUS. 
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TABLE VI: RUN SCHEDULE 


A. Rigid Model Tests, Model A in Aerolab Wind Tunnel 


Dynamic 
Test Run Pressure, Pitch Yaw 
Number (psf) Angle Angle Comments 
1 3.85 to 57.8 0 0 Ce Cw XS Noe 
in increments 
of 4,0 
Z i -30 0 : 
3 . “50 0 * 
4 if 0 30 4 
5 y - 30 30 =e 
6 " Too 30 " 
7 m 0 60 - 
8 - -30 60 : 
9 +30 60 e 
10 i 0 90 " 
mL ‘s -30 90 " 
2 me +30 90 . 
3 = 0 120 iy 
14 ‘ -30 120 ie 
15 . 50 120 "i 
16 i -30 150 i 
i? ew 0 159 : 
18 i 3 0 dis) 9, a 
19 * 0 180 oi 
20 -30 180 : 


43 


Test Run 
Number 


7h 


ee 


20 


39 


4.0 


Dynamic 


Pressure, PEeen Yaw 

(psf) Angle Angle 

Seoouto li so +30 SOs 
in increments 

Of. 0) 

a 0 210 

" -30 210 

. +30 210 

= -30 240 

iy 0 240 

ss +30 240 

m 6) 270 

- -30 270 

” +30 270 

m 0 300 

" -30 300 

" +30 300 

” 0 330 

_ -30 330 

m +30 B36 

“6 8 -30 to +30 in 0 

5 increments 

3, OD GEG Oo +90 0 

in increments 
Of 40 
46.1 16.0) Tere) ae ALG) ain 0 
5 increments 

3. SS ahem? «Sc +60 180 

in increments 
Ofearo 


Ady 


Comments 


TEN VS Np 


Cy Cu VS x 


5S, Cu VS NV = 


Cy CL Cm VS S 


CCL Cu VS Nee 


Test Run 
Number 


41 


42 


43 


44 


45 


46 


47 


48 


49 


50 


on 


3 


6 


54 


a0 


56 


Bye 


58 


3S), 


Dynamic 
Pressure, 
(psf) 


3. 835 bem O 


in increments 
Oterce( 


34862 tom. Gel 
SRG _ toea?. 3 


in increments 
of 4e0 


SU Gis) Geko Sis Gis 
3. G5mroOea2 3 
3.02mmeO 38.5 
Sai iee, Zee 


in increments 
oLaa.0 


46.1 


by eoiy cae oneal 
in increments 
of 4.0 


Pitch 
Angle 


+60 


+60 
+60 
+90 


+60 


+60 
+90 
+60 
+90 


+60 


+60 
+60 
+60 


+60 


-120 to -60 


Lia 


Lncrements 


-60 


-90 


-60 


-60 


45 


240 


210 


Comments 


“oO Cy VSN, 


Test Run 


Number 


60 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


7 


Dynamic 
Pressure 
(psf) 


46.1 


46 


Yaw 
Angle 


345mtorO 15. in 
5 increments 


O20 tom045 in 
5 increments 


OS Osram > 111 
5 increments 


OSCR to “LOS im 
5 increments 


LiOmEe 135 in 
5 increments 


140sto 165 in 
5 increments 


1ZOmGo 19521 
5 increments 


20080 8225° in 
5 increments 


2FIOVCOhS55 =n 
5 increments 


200 BECE 6) 2 in 
5 increments 


2908EOo 315 in 
5 increments 


320 to 340 in 
5 increments 


Comments 


Co Ce sae VS Y 


B. Dynamic Model Tests, West Coast Research Wind Tunnel 


Test Run Model 
No. Used 


Not Numbered 


(Note 1) B 
M-1 - M-15 

(Note 2) C 

1 - 

(Note 3) C 
mm, 11 C 

12 C 
13, 14 C 
f>, 16 C 
hi LS C 
19, 2 C 
Za 22 C 

23 C 
oy, 25 C 

24A 

(Note 4) C 

25A C 

26 C 

27 C 


Dynamic 
Pressure, 


 . Sate 
1.0' psf 


5.0 


Photographs 
(Motion Picture) 


none 


yes 


ho 


no 


yes 
film test 
sequency 16 


yes 
film test 
sequence 17 


yes 
film test 


sequence 18 


yes 


Time-History 
Trace 


none 


none 


yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
yes 
(note 5) 


yes 


yes 


yes 


Test Run 
No. 


28 


Zo 


30 


a 


a2 


Model 
Used 


C 


Dynamic 
Pressure, 


D8. 


Toh 


70 


S20 
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Photographs 


(Motion Picture) 


yes 
film test 
sequence 19 


yes 
film test 
Sequence 20 


yes 
film test 
Sequence 21 


yes 
film test 
sequence 22 


yes 
“film:-test 


csequence” 23 


yes 


no 


no 


no 


nlo 


no 


nlo 


no 


no 


no 


no 


no 


Time-History 
Trace 


yes 


yes 


yes 


yes 


yes 


yes 


yes 
yes 
yes 


yes 


yes 


yes 


yes 
yes 
yes 


yes 


yes 


Test un Model Dynamic Photographs Time-History 


No. Used Pressure, (Motion Picture) Trace 
Li=10 5 Li=1 18 6.0 no yes 
LL-12, LL-13 
LL-14, LL-15 £E 7.0 no yes 
LL-16, LL-l17 &E 8.0 no yes 
Notes: 

1. Various runs conducted in both the Aerolab tunnel and the 
West Coast Research tunnel. Model mounted on wire in order 
to observe both the yaw mode and the pitch mode. 

2. Local tie down constraints not used; these runs correspond 
to Part One of the movie film. Film test sequences 1 
through 15; film speed, 16 frames per second. 

3. Local tie down constraints were used in this and subsequent 
runs listed. 

4. Runs 24A through 33 correspond to Part Two of the movie 
film. Film test sequences 16 through 23; film speed, 

16 frames per second. 
5. Test run numbers (24A through 33, L-l through L-18, LL-l 


through LL-17) correspond with notation used on the 


time-history traces. 
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1,9 MODEL SHOWN AT ZERO PITCH AND ZERO YAW, 
2. FORCES AND MOMENTS ARE RELATIVE TO 
WIND AXIS SYSTEM (AS ABOVE), 


FliGe <2 COORDINATE AXIS SYSTEM 
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|WEIGHT | 
| 





COMPONENT WEIGHT DIA. HEIGHT 


WEIGHT A 0.30 Ibs 1.0 096 
WEIGHT 8B 0.67 Ibs 1.0 207 
WEIGHT C 0:67 Ibs 1.0 207 
WEIGHT D 0.67 Ibs 1.0 2.07 
WEIGHT E 0.67 Ibs 1.0 207 
WEIGHT F 0.67 Ibs 1.0 204 
WEIGHT G 1.45 Ibs 20 1.12 
WEIGHT #H 1.00 Ibs 20 0.77 
WEIGHT | 0.31 Ibs 1.0 0.96 
i 0.20 Ibs 
REACTOR 
CYLINCER 0.74 Ibs 
HABITAT 
CYLINCER O.EC Ibs 
MAIN BAL. AST 

TANK 1.25 Ibs 
OBSERVATION 
SPHERE 0.05 Ibs 
FIG, 8 
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Fig. ll - Model orientation: yaw zero degrees, pitch 
zero degrees; viewed from left side of model; airflow 


left to right 





Fig. 12 - Model orientation: yaw zero degrees, pitch 


zero degrees; viewed from upstream side 
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Fig. 13 - Model orientation: yaw zero degrees, pitch 


+30 degrees; viewed from left side of model; airflow 
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Fig. 
+30 degrees; viewed from upstream side 
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Fig. 15 - Model orientation: yaw zero degrees, pitch 


-30 degrees; viewed from left side of model; airflow 
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Fig. 16 - Model orientation: yaw zero degrees, pitch 


-30 degrees; viewed from upstream side 
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Fig. 17 - Model orientation: yaw 30 degrees, pitch 





zero degrees; viewed from left side of model; airflow 


left to right 
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Fig. 18 - Model orientation: yaw 30 degrees, pitch 


zero degrees; viewed from upstream side 
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Fig. 19 - Model orientation: yaw 30 degrees, pitch 
+30 degrees; viewed from left side of model; airflow 
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Fig. 20 = Model orientation: yaw 30 degrees, pitch 


+30 degrees; viewed from upstream side 
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Fig. 21 - Model orientation: yaw 30 degrees 3p itch 
-30 degrees; viewed from left side of model; airflow 
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Fig. 22 - Model orientation: yaw 30 degrees, pitch 


-30 degrees; viewed from upstream side 
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Fig. 23 - Model orientation: yaw 60 degrees, pitch 
zero degrees; viewed from left side of model; airflow 
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Fig. 24 - Model orientation: yaw 60 degrees, pitch 


zero degrees; viewed from upstream side 
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Fig. 25 - Model orientation: yaw 60 degrees, pitch 
+30 degrees; viewed from left side of model; airflow 


left to right 
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Fig. 26 ~ Model orientation: yaw 60 degrees, pitch 


+30 degrees; viewed from upstream side 
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gee? - Model orienbarion. yaw 00 degmees pitch 
-30 degrees; viewed from left side of model; airflow 
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Fig. 28 - Model orientation: yaw 60 degrees, pitch 


~30 degrees; viewed from upstream side 
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Fig. 29 = Model orientation: yaw 90 degrees, pitch 


zero degrees; viewed from left side of model; 
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Fig. 30 - Model orientation; yaw 90 degrees, pitch 


zero degrees; viewed from upstream side 
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Fig. 31 - Model orientation: yaw 90 degrees, pitch 
+30 degrees; viewed from left side of model; airflow 


left to right 
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Fig. 32 - Model orientation: yaw 90 degrees, pitch 
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+30 degrees; viewed from upstream side 
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Fig. 33 - Model orientation: yaw 90 degrees, pitch 
-30 degrees; viewed from left side of model; airflow 


left to right 





Fig. 34 - Model orientation: yaw 90 degrees, pitch 


~30 degrees; viewed from upstream side 
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Fig. 35 ~- Model orientation: yaw zero degrees, pitch 
+60 degrees; viewed from left side of model; airflow 


left to right 





Fig. 36 - Model orientation: yaw zero degrees, pitch 


+60 degrees; viewed from upstream side 
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Fig. 37 - Model orientation: yaw zero degrees, pitch 
+90 degrees; viewed from left side of model; airflow 


left to right 
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Fig. 38 - Model orientation: yaw zero degrees, pitch 


+90 degrees; viewed from upstream side 
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Pig, v= Modelworlentatlouw yaw zero degrees, patch 
-60 degrees; viewed from left side of model; airflow 
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Fig. 40 - Model orientation: yaw zero degrees, pitch 
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-60 degrees; viewed from upstream side 





Fig. 41 - Model orientation: yaw zero degrees, pitch 


-90 degrees; viewed from left side of model; airflow 
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Fig. 42 - Model orientation: yaw zero degrees, pitch 





~90 degrees; viewed from upstream side 
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Fig. 43 - Support arrangement used for yaw variation 
in five degree increments; model at zero degrees 


pitch; viewed from side; flow left to right 





Fig. 44 - Same as 43, but view From upstream side 
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Fig. 56 - Photograph of the West Coast Research 48 in. 
x 60 in. wind tunnel, showing support beams for wire 


attachment exterior to the test section 





Fig. 5/7 - The dynamic model in the test section, 


Side view 


88 






j > 
_* 
. | 


‘ 


Fig. 58 - Photograph of the dynamic model in the 





test section, view from downstream 


we.) 





Fig, 59 = Close-up of the dynamic model 
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Fig. 60 - Photograph of dynamic model instrumentation 


hardware and bearings 
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Fig, 61 - Photograph of the modified potentiometer 


used in dynamic model instrumentation 
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